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ABSTRACT 

This paper focuses on the polarized profiles of resonance scattering lines that form in magne- 
tized disks. Optically thin lines from Keplerian planar disks are considered. Model line profiles 
are calculated for simple field topologies of axial fields (i.e., vertical to the disk plane) and toroidal 
fields (i.e., purely azimuthal). A scheme for discerning field strengths and geometries in disks is 
developed based on Stokes Q-U diagrams for the run of polarization across line profiles that are 
Doppler broadened by the disk rotation. A discussion of the Hanle effect for magnetized disks 
in which the magnetorotational instability (MRI) is operating is also presented. Given that the 
MRI has a tendency to mix the vector field orientation, it may be difficult to detect the disk fields 
with the longitudinal Zeeman effect, since the amplitude of the circularly polarized signal scales 
with the net magnetic flux in the direction of the observer. The Hanle effect does not suffer from 
this impediment, and so a multi-line analysis could be used to constrain field strengths in disks 
dominated by the MRI. 

Subject headings: Accretion Disks — Polarization (Hanle Effect) — Stars: Early Type — Stars: Magnetic 
Fields — Stars: Winds — Techniques: Polarimetric 



1. INTRODUCTION 

Spectropolarimetry continues to be a valuable 
technique for a broad range of astrophysical stud- 
ies (e.g., Adamson et al. 2005; Bastian 2010), in- 
cluding applications for circumstellar envelopes. 
Advances in technology and access to larger tele- 
scopes means an ever growing body of high quality 
spectropolarimetric data. It is therefore important 
that the arsenal of diagnostic methods and theo- 
retical models in different astrophysical scenarios 
keep pace. This paper represents the fifth in a sc- 
ries devoted toward developing the Hanle effect as 
tool for measuring magnetic fields in circumstellar 
media from resonance line scattering polarization. 
The observational requirements for the effects ex- 
amined in this series are ambitious: high signal- 
to- noise (S/N) data and high spectral resolving 
power. However, these demands are being met, 
as illustrated by Harrington & Kuhn (2009a) in a 
spectropolarimetric survey of circumstellar disks 
at Ha. 

There are numerous effects that can influence 
the polarization across resolved lines. A number 
of researchers have investigated the effects of line 



opacity for polarization from Thomson scattering 
(Wood, Brown, & Fox 1993; Harries 2000; Vink, 
Harries, Drew 2005; Wang & Wheeler 2008; Hole, 
Kasen, & Nordsieck 2010). Scattering by reso- 
nance lines can generate polarization similar to 
dipole scattering (e.g., Ignace 2000a). With high 
S/N and high spectral resolution, Harrington & 
Kuhn (2009b) have identified a new polarizing ef- 
fect for lines that coincides with line absorption. 
An explanation for this previously unobserved ef- 
fect in stars is discussed by Kuhn et al. (2007) 
and is attributed to the same dichroic processes 
detailed by Trujillo Bueno & Landi Degl'Innocenti 
(1997) for interpreting polarizations in some solar 
spectral lines. Generally associated with circular 
polarization, the Zeeman effect has received acute 
attention of late as a result of techniques that co- 
add many lines (Donati et al. 1997). The method 
has been used successfully in many studies (see 
the review of Donati & Landstreet 2009). In rela- 
tion to massive stars, the technique has led to the 
detection of magnetism in several stars (e.g., Do- 
nati et al. 2002, 2006a, 2006b; Neiner et al. 2003; 
Grunhut et al. 2009). In fact, there exists a large 
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collaborative effort called the MiMeS3 collabora- 
tion (e.g., Wade et al. 2009) to increase the sample 
of massive stars with well-studied magnetic fields. 

Another method that has been used produc- 
tively in studies of solar magnetic fields involves 
the Hanle effect (Hanle 1924). This is a weak Zee- 
man effect that operates when the Zeeman split- 
ting is roughly comparable to the natural line 
broadening. Monographs that deal with atomic 
physics and polarized radiative transfer provide 
excellent treatments of the Hanle effect, includ- 
ing Stcnflo (1994) and Landi Degl'Innocenti & 
Landolfi (2004). This contribution extends a se- 
ries of theoretical papers that has enlarged the 
scope of its general use with circumstellar en- 
velopes. Building on work developed in the solar 
physics community, and using expressions for res- 
onance line scattering with the Hanle effect (e.g., 
from Stenflo 1994), Ignace, Nordsieck, & Cassinclli 
(1997; Paper I) provided an introduction of the 
Hanle effect for use in studies of circumstellar en- 
velopes. Ignace, Cassinelli, & Nordsieck (1999, Pa- 
per II) considered its use in simplified models of 
disks with constant radial expansion or solid body 
rotation as pedagogic examples. Both of these pa- 
pers approximated the star as a point source of 
illumination. Ignace (2001; Paper III) then incor- 
porated the finite source depolarization factor of 
Cassinelli, Nordsieck, & Murison (1987) into the 
treatment of the source functions for the Hanle 
effect in circumstellar envelopes. Finally, Ignace, 
Nordsieck, & Cassinelli (2004, Paper IV) calcu- 
lated the Hanle effect in P Cygni wind lines using 
an approximate treatment for line optical depth 
effects. 

The focus of this fifth paper in the series is 
the Hanle effect in polarized lines from Keplerian 
disks, of relevance for accretion disks and Be star 
disks (e.g., Cranmer 2009). As in previous papers 
of this series, the Sobolev approximation (Sobolev 
1960) for optically thin scattering lines is adopted 
for exploring the run of polarization with velocity 
shift across line profiles. This paper adopts the 
notation laid out in Paper II. The results of Pa- 
per II are expanded for Keplerian rotation with 
the inclusion of finite source depolarization (ala 
Paper III) and stellar occultation. Resonance line 
scattering polarization from such disks were ex- 
plored in Ignace (2000a) for unmagnetized disks 
without consideration of the Hanle effect. 

Drawing on these previous works, methods for 
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Fig. 1. — Illustration of the relation between 
observer and stellar coordinate reference frames. 
The observer is in the direction of Z; the symme- 
try axis for the disk is Z* . A scattering point is in 
direction f. This point lies in the equatorial disk, 
hence (p — n/2. 

computing the polarization of line profiles will be 
briefly reviewed in Section [3] In Section [3J re- 
sults for the Hanle effect from disk lines are de- 
scribed. There are three applications that will be 
considered: purely axial fields (i.e., normal to the 
disk plane), purely toroidal fields (i.e., azimuthal 
in the disk plane), and field topologies that can 
arise from the magnetorotational instability, or 
MRI (Balbus & Hawley 1991). Conclusions of this 
study and observational prospects are presented in 
Section [4] Appendices detail analytic derivations 
for special cases of the Hanle effect in Keplerian 
disks. 

2. THIN SCATTERING LINES IN KEP- 
LERIAN DISKS 

2.1. Coordinate Systems of the Model 

The focus of this paper is thin scattering lines 
from a planar circumstellar disk in which the gas 
obeys circular Keplerian motion. To describe this 
geometry and the line scattering polarization, a 
number of coordinates will be needed, which are 
introduced here. 

• A Cartesian coordinate system for the star 
is assigned to be (A*,Y*,Z*). The Z* axis 
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is the rotation axis of the star and disk. 

• Another one for the observer is (X,Y,Z). 
The observer is located at large positive dis- 
tance along the Z axis. 

• Spherical coordinates in the star frame are 
(r, i?, cp). Cylindrical coordinates are identi- 
fied by (tu, <p,Z*). 

• Cylindrical coordinates in the observer 
frame are (p, ip, Z). 

• Spherical angular coordinates defined in a 
frame of the local magnetic field at any point 
will be (6, <j>). The latter system is needed to 
evaluate the Hanle effect. 

• The Z axis is taken to be inclined to the Z* 
axis by an angle i. 

• The circumstellar disk assumed to be ax- 
isymmetric and exists entirely in the equa- 
torial plane of the star, and so it is located 
at $ = 7r/2. The scattering angle between a 
point in the disk to the observer is signified 
by X- 

Figure Q] shows a spherical triangle used in re- 
lating the observer and star coordinate systems. 

2.2. Line Velocity Shifts 

The Sobolev approximation is employed to de- 
scribe the polarimetric line profiles. The Sobolev 
approach relies on identifying "isovelocity zones" , 
which are the locus of points that share the same 
Dopplcr shift for a distant observer. The Doppler 
shift in frequency of any point in the scattering 
volume is determined by 

Az/ Z = -ful = -T— ■ (1) 

c Aui 

where v u \ is the frequency of the line transition, 
and vz is the line-of-sight velocity shift given by 

v z = -v-Z. (2) 

A Keplerian disk follows a velocity profile of the 
form 




where Rq is the innermost radius of the disk with 
tangential speed vq at that location. 
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Fig. 2. — Illustration of isovelocity zones for a Ke- 
plerian disk. The hatched region at center is the 
star, assumed to rotate ccw, as indicated by the 
arrow at top. In this example an observer is lo- 
cated in the equatorial plane. Red curves on the 
right are for redshifted velocities (dashed), with 
the smallest loop having the greatest speed. Left 
are for blueshifts (dotted). Emission at line center 
arises from the vertical solid line. 

It is convenient to define dimensionless vari- 
ables for distances and velocities. For the nor- 
malized cylindrical radius in the planar disk, w = 
r/Ro is used. The line-of-sight velocity shift then 
becomes 

vz = vqU!^ 1 ^ 2 sini sin ip. (4) 

A normalized velocity is also introduced with w z = 
v z /{v sini). 

The solution for the isovelocity zones in terms 
of w(w z , ip) is thus given by 

• 2 • 2 

sin tp sin <p 

W= — = ~ ■ ( 5 ) 

< sm z <p 

This describes a loop path that starts at Rq at 
an azimuth tpo on the front side of the disk, where 
sin tpo = w z , and then ends again at Rq and ir — po. 
The loop extends to a maximum distance at az- 
imuth p — tt/2 with a value of w = l/w%. The 
system is left-right mirror symmetric, with red- 
shifted velocities on one side and blueshifted ones 
on the other. The convention here is that the disk 
rotates counterclockwise as seen from above (i.e., 
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"prograde" ) . An example is displayed in Figure [2] 
With the arrow toward the observer as the ref- 
erence for ip, redshifted velocities come from the 
interval < <p < n, and blueshifted ones come 
from 7r < ip < 2tt. 

2.3. Line Flux 

For optically thin scattering, the observed flux 
of scattered light at a given velocity shift in the 
line is determined by a volume integral along the 
associated isovelocity loop. To describe the po- 
larization, a Stokes vector prescription is adopted, 
with standard /, Q, U, V notation. Referring back 
to Figure [TJ an edge-on disk with i = 90° would 
yield intensities with Q > and U = 0. In this 
case "North" is chosen to be along the direction 
of and Q > corresponds to oscillations of 

the electric vectors of the light being preferentially 
parallel to North. 

For unresolved disks the corresponding ob- 
served fluxes are T v = (Ti, Fq, Tjj, 0), where it 
is assumed that the Stokes-V flux is zero for the 
problem at hand. 

Generally following the notation of Paper II and 
using results from Ignace (2000a) for Keplerian 
disks, the disk is assumed to have a power-law 
surface number density given by 



S(ju) = Sq 



(6) 



for power-law exponent m and inner value Eo at 
zd = 1. Note that in the case of the Be stars, val- 
ues of m range from 3-4 (e.g., Waters 1986; Lee, 
Osaki, & Saio 1991; Porter 1999; Jones, Sigut, 
& Porter 2008), and m — 3.5 will be adopted in 
example cases. Implicit is that the lines of in- 
terest form from the dominant ion. Temperature 
and ionization variations or different disk densities 
could be included in the model, but the intent of 
this contribution is to highlight the line polariza- 
tion and Hanle effect, and so a simple power-law 
density is adopted as a baseline for the analysis. 

The Stokes fluxes for thin lines from this disk 
are given by 
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Fig. 3. — Line profiles for a disk with no magnetic 
field. Upper left: Area normalized Stokes-I line 
profiles for viewing inclinations of i = 0° (red), 
i = 30° (green), i = 45° (dark blue), i = 60° 
(light blue), and i = 90° (magenta). Lower left: 
Stokes-Q profiles plotted as the ratio of the Q- 
fiux to the scattered I-fiux as fractional polariza- 
tion (not percent). The yellow dashed line is a 
guide for zero polarization. Lower right: Like the 
Q-Profiles, Stokes-U profiles shown as fractional 
polarizations, which arise solely from the effect of 
stellar occultation. Upper right: A Q-U diagram 
across the profile. 

(hj , Hq , hjj , hy) ■ It is this factor that incorpo- 
rates the dipole scattering of resonance line po- 
larization, the impact of magnetic fields via the 
Hanle effect, and the finite depolarization factor. 
Finally, the scaling factors outside the integral are 
from equations (24) and (25) of Paper II, with 



TO 



c ; S A„i 
4Trv TO t sini' 



(8) 



where oi is the frequency integrated line cross- 
section, and 



dvj 



w v (w) |cos(/j| 



(7) 



The denominator in the integrand arises from 
the Sobolev approximation, with w v = w~ x / 2 . 
The factor | cos<p| will more conveniently be writ- 
ten as \/l — w\w in what follows. The vector 
h represents the "scattering factor" with h = 



Jo 



And 21 



(9) 



where L u ^ is the specific luminosity of the star 
and d is the distance to the star from Earth. 
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2.4. Solution for Non-Magnetic Scattering 
Polarization 

Before progressing to a consideration of the 
Hanlc effect, there is value in formulating the solu- 
tion for the line scattering polarization in the zero 
field case. This is because the Hanle effect primar- 
ily modifies the scattering polarization, meaning 
that a good understanding of the non-magnetic 
case is a necessary reference for understanding the 
Hanlc effect. Many papers have dealt with line 
formation from Keplerian disks. This work draws 
in particular on the work of Huang (1961) and Ry- 
bicki & Hummer (1983) for planar disks, and also 
fgnace (2000a) for scattering polarization without 
the Hanle effect. 

There are two main cases that are explored 
here: the simplistic case of illumination by a point 
source and then the more realistic case involving 
the effects of the finite size of the illuminating 
star. Analytic and semi-analytic solutions are de- 
scribed in the Appendices. The value of the point 
source approximation is that it gives context for 
how finite stellar size effects modify the line scat- 
tering polarization. The only portion of the point 
source approach that is reproduced here from Ap- 
pendix [X] is the expression for isotropic scattering 
of starlight by a planar disk in the Sobolev ap- 
proximation. 

For isotropic scattering of light from a point 
star, one has h = (1,0,0,0). Only the Stokes-I 
flux survives in equation (JT)). It is convenient to 
introduce a change of variable for computing the 
integration of that equation, with t = zu^ 1 . Then 
the flux becomes 

^/ = r J- x2 / dt, (10) 

Jt {w*) V* — f o 

where to =w%- The appearance of the factor of 2 
arises because the integration for the front half of 
the loop (from ifQ to ip = tt/2) is the same as for 
the back half. 

Appendix A details analytic cases for the above 
case. The addition of polarization and finite 
source effects amount to inserting multiplicative 
"weighting" functions that modify the above inte- 
grand that describes the emission contribution as 
a function of location along a loop. 

There are several finite source effects that can 
be included, such as finite star depolarization 
(Cassinelli et at 1987), stellar occultation effects 
(e.g., Fox & Brown 1991), limb darkening effects 



(Brown, Carlaw, & Cassinelli 1989), and absorp- 
tion of starlight by the disk (Ignace 2000a). Only 
two of these effects will be considered here: fi- 
nite star depolarization and stellar occultation of 
the disk. Limb darkening could be included; how- 
ever, limb darkening mainly "softens" the finite 
star depolarization factor, making it slightly less 
severe. Thus the inclusion of limb darkening does 
not seem to be a pressing issue for illustrating the 
Hanle effect in disks. 

Absorption by the disk does affect the shape of 
the Stokes-I profile; however, in the current treat- 
ment it does not influence h which determines 
the line scattering polarization. Consequently, 
absorption affects the continuum level of direct 
starlight at the line, but not the scattered Stokes 
fluxes J 7 /, J-q, or J-jj. On the other hand, a photo- 
spheric absorption line would alter the line shapes 
of the scattered flux profiles; however, results pre- 
sented here will be in terms of ratios of scattered 
fluxes, q s v = TqITj and u% = TujTi, for which 
the influence of a photospheric line will cancel. 

These ratios q% and u s v should be considered as 
polarimetric "efficiencies" . They are not generally 
what an observer would actually measure, because 
the denominator involves only the scattered light 
of the Stokes-I flux. Instead an observer would 
normally measure Tq , Tjj , and J-tot = + Tj, 
the latter being the sum of the direct starlight 
(first term) and the scattered light (second term). 
In this case q tot = Fq/{Fo + Fi) ~ Fq/Fo, and 
likewise for the Stokes-U polarization, assuming 
that the scattered flux in the line is small com- 
pared to the direct flux by the star itself. As a 
result, (q to t, "tot) = (ql,ul) x (Fi/Fo)- It is worth 
noting that the q s v and u s v "efficiency" profiles are 
independent of the line optical depth To in the thin 
limit, since J 7 /, Fq, and Tjj scale linearly with To- 

Equation IB II gives the scattering function com- 
ponents of h in the case of no magnetic field with 
the star treated as a point star of illumination. 
Modifications to those functions that allow for 
scattering by a finite sized and uniformly bright 
central star were determined in Paper III. For the 
case of zero magnetic field, the scattering function 
is shown explicitly with factors arising from the 
finite stellar size. These factors also appear in the 
scattering function with the Hanle effect in the 
same positions as without the Hanle effect. The 
vector components of h are: 
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h= < 



r hj 

hu 



W + ^E 1 ^ [A (1 -3 cos 2 i) 
+ sin 2 i cos 2ip\ 
| -Ei /i» [sin 2 i 

cos 2 cos 2<£>] 



1 



\Ei |U» cosi sin 2y>, 



(11) 



where 



W> = 4w 2 PF(ot), (12) 
with the dilution factor W(ru) = 0.5 (1 — /i*), and 



(13) 



is the finite depolarization factor when there is no 
limb darkening. The factors W and /i» represent 
corrections to the scattering function that account 
for the effects of finite star size in relation to the in- 
cident radiation field at a scattering point around 
the star. 

Scattering by resonance lines can be approxi- 
mated as part dipole and part isotropic. The pa- 
rameter Ei is a factor with values from to 1 that 
represents the extent to which a resonance line 
scatters like a dipole radiator (see Chandrasekhar 
1960). Only the dipole portion of the scattered 
light contributes to observable polarization. In the 
solar literature, it is more common to use the no- 
tation W2 for the fraction of scattered light that 
is dipole- like (e.g., Stenflo 1978). 

Accounting for stellar occultation breaks the 
back-front symmetry of the isovelocity loops. The 
front loop still has lower and upper limits of to and 
1, but the back half has limits of to and t* < 1. 
Here t* is associated with the minimum value of 
w = w* at the back-projected limb of the star 
corresponding to a maximum azimuth of (p = tp* . 
An expression for is given in Ignace (2000a), 
expressed here as 



t* = \J 1 — sin 2 i cos 2 ip* = \J cos 2 i — w% sin 2 i /t*. 

(14) 

This expression provides a cubic relation in i* as a 
function of fixed viewing inclination i and velocity 
shift w z . 

Numerical results for polarized profiles without 
the Hanle effect are shown in Figure [3] This Fig- 
ure and the ones that will follow use m = 3.5 and 
Ei = 0.5. Figure [3] should be compared to Fig- 
ure [9] for the point star case. For point star illu- 
mination, the maximum polarization occurs at the 



line wings and J-jj = by symmetry. This changes 
dramatically when the effects of the finite stellar 
size are included. 

In Figure 131 each panel has 5 curves color-coded 
for different viewing inclinations: sin 2 i — is red, 
sin 2 i = 0.25 is green, sin 2 i = 0.5 is deep blue, 
sin 2 i = 0.75 is light blue, and sin 2 i = 1.0 is ma- 
genta. The values correspond to viewing inclina- 
tion angles of i = 0°, 30°, 45°, 60°, and 90°, re- 
spectively. The upper left panel displays Tj pro- 
files that are normalized to have unit area. In the 
lower right panel for the q s v polarization, the fi- 
nite depolarization factor shifts peak polarization 
to lower velocity shifts as compared to the point 
star case, and q s v remains symmetric about line 
center. Stellar occultation leads to the survival 
of a small u s v polarization as shown in the lower 
right panel. The upper right panel is a q%-u s u plot 
across the respective line profiles. Note that the 
scales in q s u and u s v differ by about an order of mag- 
nitude. These small loops signify position angle 
rotations across the observed polarized line pro- 
files, with polarization position angle tpp given by 
tan2?/>p = ul/ql- 

We next turn our attention to the Hanle effect. 
The results of this section will prove valuable in 
following how the Hanle effect modifies the run of 
polarization across the line profile. 

3. THE HANLE EFFECT IN KEPLE- 
RIAN DISKS 

The Hanle effect applies to resonance line scat- 
tering and can be interpreted in semi-classical 
terms as a precession of an oscillating emitter that 
occurs over the radiative lifetime of the line emis- 
sion. An angular quantity can be defined to rep- 
resent the effective amount of precession, as given 

by 



tan ct2 



2 9l u l 



B 



Bn 3 



(15) 



where ujl is the angular Larmor frequency, A n \ is 
the Einstein radiative rate for the transition of in- 
terest, <?l is the Lande factor of the upper level, B 
is the magnetic field strength in Gauss, and -BHan 
is the Hanle field sensitivity defined as 



£?Han = 56.9 — Gauss, 



(16) 



with j4g the radiative rate normalized to 10 9 Hz. 
The Hanle effect manifests itself in the scatter- 
ing physics with the appearance of cosine and sine 
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functions of the angle oli (and a related quan- 
tity ol\ = 0.5a2) within the scattering functions 
h. Note that with no magnetic field, B = and 
q?2 = 0. The limit of a high Larmor frequency 
yields 012 ~ tt/2. In this latter case, the Hanlc 
effect is said to be "saturated" . 

There are a few rules of thumb to help under- 
stand the Hanle effect in polarized lines. First, 
modifications to the polarized profile will depend 
on the orientation of the field in relation to both 
the direction of incoming radiation and outgoing 
radiation (i.e., the observer). Second there is no 
Hanle effect when the incoming radiation field is 
symmetric about the field direction. These first 
two "rules" indicate that there is no Hanle effect 
for radial magnetic fields (assuming a spherically 
symmetric star and no diffuse radiation) and that 
only non-radial components of the field contribute 
to modifying the polarization. However, the third 
point is that in spite of this, the angle a 2 is sen- 
sitive to the total field strength: both radial and 
non-radial components, at the site of scattering. 
Fourth, and finally, the Hanle effect is sensitive 
only to the field topology in the saturated limit, 
not the field strength. 

With these rules of thumb, consideration of the 
Hanle effect in three particular cases follow: a 
purely axial field, a purely toroidal field, and last 
a scenario involving fields as arising from the op- 
eration of MRI in disks. 

3.1. Axial Magnetic Field 

In this case an axial magnetic field is envisioned 
as threading the disk perpendicular to the equato- 
rial plane, hence B — Bz, (n?) Z*. As for the field 
distribution throughout the disk, the following is 
adopted as a conservative limiting case: 

B Zt =B vj-\ (17) 
and so the Hanle angle 012 is determined by 

tana2 = t; nj" 1 = bot. (18) 

where b — -Bo/£?Han for B a field strength at 
the inner disk radius. It is important to note that 
if Bo <C -Brian, then the Hanle effect is weak ev- 
erywhere in the disk. If B n >• £>Han, then the 
inner disk will be in the saturated limit; however, 
there will be a transition to a weak Hanle effect 
at some radius in disk, characteristically where 
w = Bo/Bnan- As a general rule, for a given 
magnetized disk, different lines will have different 




-1 -0-5 0.5 1 -1 -0.5 0.5 

Velocity Shift, w. Velocity Shift, w. 



Fig. 4. — The Hanle effect for a disk threaded by a 
purely axial magnetic field. Left is for q s y and right 
is for u s v . Profiles are shown for different viewing 
inclinations with i = 30° at top, then i = 45°, 
i = 60°, and i = 90° at bottom. The colors are for 
different field strengths with = 0.01 (red), 0.1 
(green), 1.0 (dark blue), 10 (light blue), and 100 
(magenta). 

values of i?Han, and thus will be sensitive to dif- 
ferent aspects of the field at different locations in 
the disk. 

Expressions for h in the point star limit are 
given in Appendix [CJ Revisions to those expres- 
sions for finite star depolarization and stellar oc- 
cultation is the same as in equation (fTTj). Results 
are shown in Figure 2J Profiles of q s v and u% are 
displayed at left and right, respectively. Panels 
from top to bottom are now for different view- 
ing inclinations of sin 2 i — 0.25,0.5,0.75, and 1.0. 
Different colors are for different values of bo with 
b Q = 0.01 (red), 0.1 (green), 1.0 (dark blue), 10 
(light blue), and 100 (violet). 

The Stokcs-Q profiles are symmetric about line 
center, whereas the U profiles evolve from weak 
and anti-symmetric, owing to stellar occultation, 
to strong and mostly symmetric, to a null profile 
at the strongly saturated limit. In the point star 
case, the saturated limit yields u s u = and q s v = 
constant for all velocity shifts. With finite source 
effects, qf, is mostly flat throughout the central 
portion of the profile, but goes to zero in the wings. 
In the presence of limb darkening, the q s u profile 
would become more flattened toward the wings. 
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Fig. 5. — A Q-U diagram for the polarized profiles 
shown in Figure 0] with the same color scheme in 
relation to results for different values of bo- The 
strongly "saturated" case of b = 100 is present in 
magneta. This case has a modest range in q s v from 
to about 0.1, but with u% everywhere, the 
profile is only a short horizontal line in the Q-U 
plane and difficult to see in this figure. For the 
edge-on case, u s v = and all profiles degenerate to 
horizontal lines from to a maximum value of q„ 
that depends on bo- 

Figure [5] displays q s v -u% curves across the line 
profile. Note that the two axes have different 
scales: variations in u s v are smaller than for q s y . 
The colored curves are the same cases as shown in 
Figure |4j Viewing inclinations are indicated. For 
the edge-on case of i = 90°, the curves are all de- 
generate with u s v — 0. The red loop is for a very 
weak field and is seen to be top-bottom symmet- 
ric in this space. As the field is increased, q s u and 
u s v profiles become individually more nearly sym- 
metric. In the ql,-u a v space, this results in curves 
with only small loops, for example the modest field 
case with the dark blue curve. As the disk enters 
the saturated limit, u s v drops to zero, q s u becomes 
somewhat flat-topped in appearance (again, ex- 
cept at the line wings), which degenerates mainly 
to a point in the q s v -u% space for most velocity 
shifts, with extension to zero polarization only for 
the line wings. 



3.2. Toroidal Magnetic Field 

A completely azimuthal field configuration has 
also been considered. As in the axial field case, the 
toroidal field strength is assumed to decrease in- 
versely proportional to w, and so equation (|T8|) 
remains valid. This means that differences in 
the model polarized profiles between the axial 
and toroidal field configurations arise because the 
Hanle effect is sensitive to the vector field orienta- 
tion. 

As with previous considerations, some analytic 
results are derivable in the point star limit, and 
these are detailed in Appendix |TJ] for a disk with 
a toroidal field. With a toroidal field, the geome- 
try associated with determining the Hanle effect is 
more complex than in axial field case. Geometrical 
relationships between the various angles defining 
the scattering problem with a toroidal field in a 
disk were given in the Appendix of Paper II and 
will not be repeated here. Results for the calcu- 
lation of polarized profiles, including the effects 
of the star's finite extent, are shown in Figure [6] 
This Figure is presented in the same manner as 
Figure H] for the axial field case: q s v is at left and 
u s v at right; the panels are for different values of 
sin 2 i; the colors are for different field strengths 
characterized by bo- 

Note the marked differences in the line polar- 
ization between the toroidal case and that of an 
axial field. The Stokes-U profiles are always anti- 
symmetric. Like the axial case, q s y profiles are sym- 
metric, but the behavior is quite different. For ex- 
ample the edge-on disk case at lower left indicates 
that the profile polarization completely changes 
sign, from positive definite at every velocity shift 
when B = to everywhere negative in the satu- 
rated limit. This implies a 90° rotation in position 
angle between these two limiting cases for the en- 
tire profile. For intermediate field strengths, there 
is a position angle rotation across the profile that 
occurs at different velocity shifts. In the axial field 
case, there is never a position angle rotation for an 
edge-on disk. For the q%-u s v shapes, differences in 
the polarizations are especially clear, with results 
for the toroidal field in Figure [7] to be compared 
against those for an axial field in Figure [5] 

What is the source of these differences between 
the axial field and the toroidal one? The field ori- 
entation with respect to the viewer is clearly the 
key to the interpretation. For the axial case, the 
field at every point in the disk has some compo- 
nent of the magnetic vector directed toward the 
observer (if seen at i < 90°) or away from the ob- 
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Fig. 6. — Similar to Figure[4]but now for a toroidal 
field instead of an axial one. From top to bottom, 
panels are for the same viewing inclinations, and 
the color scheme is for the same values of bo. The 
ql profiles are symmetric although differing in de- 
tails from the axial field case. More striking is 
that the u s v profile is decidedly antisymmetric for 
the B v case, whereas profiles tend toward being 
symmetric in the axial field case. 

server (if seen at i > 90°). Thus, an axial field 
leads to a net projected magnetic flux of one sign 
or the other, which is true for every isovelocity 
loop. The toroidal field is manifestly different, 
since one side of the disk has components toward 
the observer and the other side has them away. 
For an axisymmetric toroidal field, the projected 
net magnetic flux is identically zero for the entire 
disk, but is oppositely signed in isovelocity loops 
for blueshifted velocities versus redshifted ones. In 
terms of the semi-classical precession description, 
flipping the field by 180° amounts to a preces- 
sion in the opposite direction. In the axial field 
case, the precession of the radiating oscillator is 
uniform - entirely clockwise (cw) or counterclock- 
wise (ccw). Both cw and ccw precessions occur for 
a toroidal field configuration, with one precession 
occurring in half the line profile, and the opposite 
for the other half. 

To illustrate this effect, Figure [8] shows line pro- 
file results for a toroidal field that now goes in the 
other direction as compared to Figure [6l The ql 
profiles are nearly the same, but the u s v profiles are 
nearly mirror images of one another, slight differ- 



Fig. 7. — Similar to Figure[5]but now for a toroidal 
field instead of an axial one. The combination 
of antisymmetric u s v with symmetric q% leads to 
Q-U loops that are top-bottom symmetric in this 
space. Note that the axis scale is different for u s u 
as compared to q s u . 

ences are a result of the stellar occultation, which 
does not flip when the field orientation is reversed. 

3.3. Magneto-Rotational Instability 

A large literature has developed over the last 
twenty years in relation to the magnetorotational 
instability (MRI; Balbus & Hawley 1991, 1998). 
There is an interesting history about this effect 
(e.g., Balbus 2003). The instability can be il- 
lustrated through an analogy to two masses in a 
differentially rotating disk that are slightly offset 
from each other along a radius. These two masses 
are connected by a spring, to represent the effect 
of an axial magnetic field. The end result is that 
the two masses on different orbits seek to increase 
their displacement from one another. Coupling by 
the spring leads to a runaway situation ensues. 

Key for the context of magnetic diagnostics is 
how this instability impacts the field topology and 
strength throughout the disk. The mass-spring 
analogy above has built into it the assumption 
of flux freezing. The separation of the masses 
along with the differential rotation would appear 
to evolve the axial field through the disk into 
a toroidal one. Different researchers have stud- 
ied the operation of the MRI in accretion disks 
through both semi-analytic work and numerical 



9 




-1 -0-5 0.5 1 -1 -0.5 0.5 

Velocity Shift, w. Velocity Shift, w. 



Fig. 8. — The same cases shown in Figure [6] except 
now for a toroidal held with the opposite sense of 
rotation. Note how the u s v profiles are essentially 
the reverse of those shown in Figure [5] 

simulations (e.g., Balbus & Hawley 1991, 1992; 
Hawley & Balbus 1991, 1992; Hawley, Gammie, 
& Balbus 1995; Hawley 2000; Fromang & Stone 
2009; Lesur & Longaretti 2009; Maheswaran & 
Cassinelli 2009). 

The goals of the MRI simulations are to un- 
derstand better the physics of angular momentum 
transport through disks and disk structure. How- 
ever, this paper seeks better insight into whether 
and how disk magnetism might be directly de- 
tected, with a focus on the Zeeman and Hanle ef- 
fects for spectral lines. Ignace & Gayley (2008) 
reported on a simplistic calculation of the Zeeman 
effect and the Hanle effect for a Kcplerian disk 
with a purely toroidal field. Here, application of 
the Hanle effect to Keplerian disks has been de- 
veloped in a more complete way. But before dis- 
cussing its application to the MRI scenario, it is 
worth commenting on the conclusions of Ignace & 
Gayley for the use of the Zeeman effect. 

The MRI leads to a field topology that consists 
of a toroidal component and a randomized compo- 
nent. The toroidal components exists in mixed po- 
larity, by which it is meant that some sectors run 
cw and others run ccw. For the weak longitudinal 
Zeeman effect, the Stokes- V flux scales with the 
net magnetic flux associated with a spatial resolu- 
tion element (e.g., Mathys 2002). With bulk mo- 
tions frequency (or wavelength or velocity) resolu- 



tion ultimately maps to geometrical zones at the 
source (an example of this involving the Sobolev 
approximation was detailed by Gayley & Ignace 
2010 for spherical winds). The key point is that if 
both the randomized field and the polarity of the 
toroidal component changes on small scales, then 
the flux of circularly polarized light arising from 
the Zeeman effect will be strongly suppressed ow- 
ing to little net magnetic flux. 

The issue of variations of the field on "small 
scales" must be evaluated with care. In the 
Sobolev approximation for purely Keplerian rota- 
tion, we have seen that the isovelocity zones are 
"loop" structures. These loops extend out to large 
radius for low velocity shifts, but they can be quite 
small at high velocity shifts. The highest velocity 
shifts degenerate to a pair of points for either limb 
of the star at the projected stellar equator! More 
importantly, every isovelocity zone intersects with 
the innermost radius of the disk. (If the disk ex- 
tended down to the star, this would be the photo- 
spheric radius.) A steep power-law density ensures 
that the bulk of emission or scattered light comes 
from inner disk radii. As a result, the Zeeman ef- 
fect would be most sensitive to a magnetic field 
at these locations, and thus "small scales" refers 
to turnovers in the field direction that are small 
compared to the inner radius of the disk. Conse- 
quently, the detection of the Zeeman effect from 
a disk with a given spectral resolution and den- 
sity distribution constrains the characteristic spa- 
tial wavelengths at which the field turns over with 
radius from the star, azimuth around the star, and 
vertical height through the disk. 

A detection of the Zeeman effect in a disk has 
previously been reported by Donati et al. (2005) 
in the case of FU Ori. This is thought to be an 
accretion disk with a disk wind, as evidenced by 
some lines showing P Cygni absorption. However, 
the circular polarization profile is argued as being 
associated with the innermost region of the rotat- 
ing disk where the field is of kilogauss strengths. 
The detection implies a net magnetic flux per spec- 
tral resolution element, and thus sets limits on the 
turnover (or "tangledness" ) length scale for the 
disk field, if indeed the MRI is operating in this 
case. 

It is interesting to consider signals that could 
result with the Hanle effect. The Hanle effect can 
operate in regions where magnetic fields are "tan- 
gled" or randomized. This means that spatial av- 
erages of (B) tend toward zero although (B 2 ) does 
not, such as is the case for the MRI mechanism. 
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An extensive literature exists for the Hanle effect 
with random fields in applications to solar studies 
(e.g., Frisch et al. 2009, and references therein). 
Here I simply want to outline some of the limiting 
behavior in applications to disks. 

Consider a Keplerian disk that contains ev- 
erywhere a truly randomized magnetic field. If 
the field is weak at all locations, meaning that 
B <C -BHan, then of course a <f v line profile results 
as in the case of no Hanle effect. But if the field 
is strong, such that the denser regions of the disk 
are largely saturated, then the behavior is much 
different. With different field orientations, one ex- 
pects that u% will yield a null profile, by symmetry 
considerations. However, the q s v profile is different. 

As a specific example, consider the resulting 
Tq line from an edge-on disk. Without a field, 
the polarization at line center would be zero, ow- 
ing to forward scattering of unpolarized starlight. 
With a randomized field, the polarization will still 
tend toward zero at line center. In the point star 
limit, some net polarization is expected to survive 
in the line wings. This polarization will be sig- 
nificantly reduced in comparison to the zero field 
case. If Ei were unity, the polarization at the line 
wings would be 100%, since the scattering geome- 
try is 90° . In analogy with considerations of scat- 
tering polarization off the solar limb, a reduction 
in polarization by a factor of 5 for isotropically 
distributed fields should be expected (see Stenflo 
1982). 

Now consider the introduction of a sustained 
toroidal field component. Of course toroidal fields 
were considered in the previous section. Now how- 
ever, the toroidal field has polarity flips within the 
disk - B v is alternately cw or ccw at essentially 
random points within the isovelocity zones. What 
this means is that there is a sign change in the 
direction of Larmor precession in the classical pic- 
ture of a harmonic oscillator. The effect of this 
is to drive the u s v signal to zero faster than if the 
toroidal field had one sense of polarity. In the sat- 
urated limit, the q s v profile is unchanged, because 
the surviving polarized signal does not depend on 
polarity at all. This is quite different from the 
Zeeman effect that is sensitive to the net magnetic 
flux. For the Zeeman effect, the circular polariza- 
tion will be suppressed when B v switches polarity 
on small spatial scales. 

4. CONCLUSIONS 

Polarized line profile shapes from magnetized 
Keplerian disks have been calculated under a num- 



ber of simplifying assumptions: the disk is geo- 
metrically thin; the scattering lines are optically 
thin; primarily simple fields were considered (axial 
or toroidal); and no account was taken of photo- 
spheric absorption lines. On the other hand, the 
model profiles do include finite source depolariza- 
tion and the effects of stellar occultation. The 
presentation of results focused on the polarimet- 
ric "efficiencies" q s v and u% with a description of 
how to identify the occurrence of the Hanle effect 
in scattering lines from disks. In addition, a dis- 
cussion was presented for the Hanle effect from a 
magnetized disk in which the MRI mechanism is 
operating. 

One of the main conclusions from this work is 
that axial and toroidal fields in disks are easily dis- 
tinguishable through an analysis of q s v -u% figures 
for the run of Stokes polarizations across line pro- 
files. Although a u s u profile does exist even without 
the Hanle effect, owing to stellar occultation, its 
amplitude is quite small. The strongest u s v pro- 
files result when much of the inner disk, where 
most of the scattered light is produced, has val- 
ues of B/Bn an of order a few. If the inner disk is 
mostly in the saturated limit of the Hanle effect, 
the u s v profile becomes a null profile for both axial 
and toroidal fields. 

One should bear in mind that the polarized ef- 
ficiencies are upper limits to the polarizations that 
would actually be measured, since the efficiencies 
are with reference to the scattered light only and 
do not take account of the direct starlight from 
the system. Since this starlight is expected to be 
largely unpolarized, its contribution acts to "di- 
lute" the polarization substantially below the ef- 
ficiency levels reported here. For example, if the 
line is relatively weak at 20% of the continuum 
level, then the expected measured polarizations 
would have fractional values at about 1/5 of the 
efficiency values, resulting in line polarizations of 
around l%-2% for q s v and 0.5% or less for u s v . 

As percent polarizations, such values would ap- 
pear to be easily measurable, but in practice there 
are several challenges. First, a spectral resolu- 
tion yielding several points across the polarized 
profile is needed. Circumstellar disks have ro- 
tation speeds of order 500 km s^ 1 . A requisite 
velocity resolution of perhaps 50 km s _1 would 
then be needed, implying a resolving power of 
A/ A A w 6000. Harrington & Kuhn (2009a) have 
demonstrated that such resolving powers can be 
achieved in spectropolarimetry; however, the next 
requirement is that of finding suitable scattering 
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lines. 

The very interesting effects seen in the sam- 
ple of Harrington & Kuhn (2009a) exploits a rel- 
atively new effect of enhanced polarization that is 
coincident with regions of higher line absorption, 
a consequence of optical pumping effects. For the 
Hanlc effect, or even for non-magnetic resonance 
scattering, lines that are predominantly scatter- 
ing are needed. For hot star disks, such as the 
disks of Be stars, resonance scattering lines are 
generally to be found at UV wavelengths. This re- 
quires space-borne spectropolarimeters. Although 
the Wisconsin Ultraviolet Photo-Polarimetcr Ex- 
periment (WUPPE) obtained exciting new results 
from UV polarimetry, its resolving power was 
only about 200 (Nordsieck et al. 1994). A new 
UV spectropolarimeter called the Far-Ultraviolet 
SpectroPolarimeter (FUSP) is a sounding rocket 
payload that will have a resolving power of about 
1800 (Nordsieck 1998). Although possibly two low 
for circumstellar disks, it will be suitable for study- 
ing scattering line polarizations from high velocity 
stellar wind sources. 

As a matter of practical analysis, how should 
spectropolarimetric data best be managed to mea- 
sure a Hanle effect? Plotting the velocity shifted 
polarizations in the q s „-u s v space appears to be most 
promising. The sequencing of the analysis for scat- 
tering lines from disk sources might proceed as fol- 
lows: 

• Subtract off the continuum polarization in 
the vicinity of the line of interest. This will 
counter the effects of both interstellar polar- 
ization and any other broadband source po- 
larization, such as may arise from Thomson 
scattering in the disk. 

• Plot the Q and U line fluxes, not fractional 
or percent polarizations that would be de- 
rived through normalization by the total 
I-flux. The reason for plotting polarized 
fluxes is that many common UV resonance 
lines are doublets (e.g., Nv, Siiv, and Civ 
UV resonance doublets). The shorter wave- 
length component ("blue") has E\ = 0.5, 
and the long wavelength component ( "red" ) 
has Ex = (e.g., Tab. 1 of Paper II). If the 
lines are thin, then the polarized flux from 
the red line will not be influenced by the blue 
one. However, if the lines are sufficiently 
closely spaced, then Fi will be a blend, and 
normalization by that blend would artifi- 
cially skew the Q-U figure shape, making 



interpretation more difficult. If the doublet 
components are well separated, then relative 
polarizations could also be used in what fol- 
lows. 

• Determine whether the resultant figure for 
the line polarization shows any axis of sym- 
metry. If so, then either (a) there is no Hanle 
effect, or (b) there is a Hanle effect with a 
dominant toroidal component. For an axis 
of symmetry, a rotation of the figure from 
observer Q-U axes to a source set of axes 
Q'-U' could be accomplished with a Mueller 
rotation matrix. After rectifying the figure 
in this way, the relative amplitude of Tjj 
should be compared to Tq. If Tjj <C Tq, 
then there is likely little Hanle effect or the 
disk is in the saturated limit of the Hanlc 
effect. If Tjj ~ Tq, then a Hanle effect is 
required with 6 ~ few in the disk where 
the bulk of scattered light is produced. This 
means B rss Bn an . 

• If there is no symmetry axis to the Q-U fig- 
ure, then a Hanle effect involving an axial 
field is most likely the culprit. Recall that a 
radial field could be present, but this would 
give no Hanle effect by itself. 

• For identifying a field distribution arising 
from MRI, things are more complicated. If 
the Hanle effect is in operation, then the 
Stokes-U flux is likely driven to zero, even 
if &o ~ 1- A net Tq profile should survive; 
however, it may not be symmetric. Using 
an oversimplification, each isovelocity zone 
can be considered to have a different effec- 
tive bo value. This is already the case in 
the pure axial or toroidal field case that is 
axisymmetric, but the variation of the effec- 
tive bo with velocity shift is monotonic (in 
a flux- weighted sense). With MRI, one ex- 
pects non-monotonicity from fluctuations of 
the field strength. This amounts to the in- 
troduction of amplitude fluctuations across 
the polarized line. 

There are complications to the above approach. 
For example, in the case of the Be star disks, 
there is substantial evidence for one-armed spiral 
density wave patterns that make the disk density 
non-axisymmetric (e.g., Okazaki 1997; Papaloizou 
& Savonije 2006). For thin Thomson scattering, 
there is little observational consequence of this ef- 
fect which is mainly a redistribution of disk scat- 
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terers in a point antisymmetric pattern (e.g., Ig- Foundation (AST-0807664). 

nace 2000b). However, with resolved line pro- 
files, the non-axisymmetric density distribution 
will produce asymmetry in the q s v profile and will 
be a new source of net u% signal, also asymmet- 
ric. Certainly it will be important to model the 
Ti profile self-consistently along with the polar- 
ized profiles. 

Note that there may be some concern about 
the choice adopted for the radial distribution of 
the field strength through the disk. The choice of 
\B\ oc w^ 1 is perhaps the most shallow distribu- 
tion that one could expect. A steeper gradient of 
the field strength will restrict the operation of the 
Hanle effect to a more restricted range of radii in 
the disk. If the interval in radius where bo > 1 
becomes narrow in relation to where most of the 
scattered light is produced, then the Hanle effect 
will be irrelevant for the observed polarization. 

What are the next steps in formulating better 
diagnostics of the density and magnetic field struc- 
ture in disks? Most of work considered in this se- 
ries of papers has focused on optically thin lines in 
an attempt to gain a better understanding of how 
the Hanle effect influences line polarizations that 
form in circumstellar media. In Paper IV the is- 
sue of optical depth effects for P Cygni lines from 
stellar winds were considered by treating regions 
of high optical as contributing no polarization at 
all. Although simplistic, insight was gained into 
how optical depth effects provide an additional 
spatial "filter" in terms of where bulk of line po- 
larization will be produced. Naturally, rigorous 
radiative transfer and a more realistic disk model 
(i.e., not simply planar) is badly needed to extend 
the considerations of this paper. Even with thin 
line scattering, it would be useful to explore how 
non-axisymmetric disk models, such as the one- 
armed spiral density wave pattern for Be disks, 
will modify the line polarization. Moreover, winds 
driven off magnetized disks (e.g., Konigl 1989; 
Proga, Stone, & Drew 1998; Proga, Stone, & Kail- 
man 2000) have been ignored entirely in this work. 
These are all new calculations that will need to be 
considered in the future. 

The author is indebted to the anonymous 
referee whose comments have improved this 
manuscript. Appreciation is expressed to Ken 
Gayley, Joe Cassinelli, and Gary Henson for en- 
gaging conversations about line polarizations and 
the Zeeman and Hanle effects. This work was 
supported by a grant from the National Science 
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Fig. 9. — Analytic line profile shapes from isotropic scattering for the density exponent parameter m = 3 
at left and m = 4 at right. The flux profiles are shown at top and qf, profiles at bottom. The different 
curves are for different viewing inclinations, with sin 2 i — 0.0, 0.25, 0.5, 0.75 and 1.0. In each case E\ = 0.5. 
The flux profiles are normalized to the total line flux in the isotropic case. The normalized isotropic line is 
plotted in light blue in the upper panels. 
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APPENDIX: SPECIAL CASES FOR THE LINE PROFILE SHAPES 

In these Appendices, a number of instructional cases are considered for the polarized line profiles shapes 
from a Keplerian disk when the illuminating star is treated as a point source. This means that both stellar 
occultation and the finite star depolarization factor are ignored. A consequence of this approximation is that 
a non-zero u s v profile can only result from the Hanle effect. Before considering polarized line profiles, Stokes- 
I profile shapes are derived for the case of isotropic scattering. These solutions form the base emissivity 
function from which the polarized lines are constructed. 



A. The Case of Isotropic Scattering 

Isotropic scattering corresponds to E± =0, and it means there is no polarization from resonance line 
scattering. Of course, that also means there is no Hanle effect, regardless of the field strength. 

Even though there is no Hanle effect, the isotropic case is useful to explore as a reference for the production 
of the Stokes-I line shape. The integrand for the line emission as a function of velocity shift in the observed 
line represents the contributions by the disk density and the Sobolev effect for the profile shape. Allowing 
for Ei ^ and the Hanle effect simply represents new weighting functions for non-isotropic scattering that 
multiply the integrand from the isotropic case. 

The flux of line emission at normalized Doppler shift w z is 

/■l j.m— 1 

H so = T E a x2 ~j==dt, (Al) 

where the factor of 2 arises from the back-front symmetry of the integration along the isovelocity zone. As 
a reminder, t — zu^ 1 and to = 

Again the preceding expression is only valid in the point star approximation. The power law exponent m 
is from the surface density distribution that is assumed to be a power law of the form £ = E w~ m . This 
formulation leads to symmetric double-peaked line profile shapes for m > 2. Larger values of m result in 
profiles that have more pronounced double-horns at greater velocity shifts from line center. 
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For integer values of to, the integral is analytic, and solutions for to — 3 and m — 4 are given here by way 
of example. For m = 3 the result is 



J™ (3) = t J x 2 / -=£= dt, = r Jo x A (3 + 4u£ + 8w z 4 ) a/i-w 2 , (A2) 

and for m = 4, 

J*°(4) = r„J„x2| -^== dt, - r Jo x 1 ^/l-to* + | ™ z 2 J*°(3). (A3) 

Note that increasing values of m also result in profile shapes of lower amplitude. In fact it is possible to 
solve for the integrated light from a resonance scattering line in some special cases. Defining 

r +1 ■ 

F tot (m) = J ^ F™dw z , (A4) 

a few selected results are F tot (l) = 47r • t J , F tot (2) = 5ir/3 ■ t Jo, F tot (3) = 2ir/3 ■ r J , and F tot (4) = 
ir/2 ■ r Jo- 
lt would appear that with to = constant, different values of F% %, and thus lines of different brightness 
levels, could result. This seems counterintuitive if lines characterized by different values of to have the same 
optical depth. However, tq depends only on the density at the inner portion of the disk Soj not the value of 
to. Thus when the line is optically thin, the appropriate optical depth to use is one that is angle averaged, 
similar in spirit to f in Brown & McLean (1977) for optically thin electron scattering. Hence use of a new 
optical depth parameter T; — To x tq J 7 a/F to t (to) would ensure that lines of different to values with isotropic 
scattering will have the same total line emission (i.e., "area under the curve") even though they have different 
profile shapes. 



B. The Case of Resonance Scattering with B = 

For resonance line scattering with E\ ^ but with B = everywhere, the vector scattering function h 
greatly simplifies. Following Paper II, we have that 6 = ip, a-2 = 0, C = cos 2cp, D = sin2y>, and tf) s = 0. The 
components of the phase function become 

!hi = 1 + §.Ei [| (l - 3cos 2 i) +sin 2 i cos2<^] 
h Q = |£a [sin 2 i- (l + cos 2 i) cos 2tp] (Bl) 
hjj = g-Ej cosz sin2<ys. 

The Stokes flux of scattered light is 

/■l j.m — 1 

P c = r Jo x 2 / -== halt. (B2) 

Jt V t - h 

Note that care must be taken in dealing with terms that are odd and even in if between angles of and tt. 
For example, hjj is odd around the loop, ensuring that Txj = 0. Accounting for the odd/even terms, and 
using the fact that sin 2 if = t /t = w%/t, solutions for the scattered flux in Stokes-I and Stokes-Q are: 



J/ sc = (1 + \E X - lE lC os2i) Ji so (m) - \E X sin 2 i t Pj°(m - 1), (B3) 
T Q SC = -lEi cos 2 i Ji so (TO) + |Fi (1 + cos 2 i) t Ji so (m - 1). (B4) 

Note at the extrema of the line wings, J^ so (to) = J^ so (m') and to = 1, and q s v = 3i?i/(4 — E\) always. 

The resultant polarization across the line is entirely in Stokes-Q when the observer's reference axes are 
aligned with the symmetry axis of the star. With E\ = 0.5, Figure [9] shows a plot of J| c for to = 3 (left) and 
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m = 4 (right) at different viewing inclination angles of i — 0°, 30°, 45°, 60°, and 90° along the top panels. 
The profiles are normalized with respect to the total emission produced if the line had been isotropically 
scattering. At bottom is the relative fractional polarization qf, — Fq/J-f^ for m = 3 and m = 4 for the same 
viewing inclinations. Note that at the edges of the line, q s v = 3/7 for E\ = 0.5, independent of the viewing 
inclination, as expected under the point star approximation. 

The goal here is to illustrate the polarimetric "efficiency" . The actual measured fractional polarization 
would be much smaller owing to dilution by direct starlight. These efficiencty curves are relatively smooth 
functions of velocity shift. This smoothness is partly due to the fact that dipole scattering is a fairly slowly 
varying function of location around the disk and also because isovelocity loops sample a range of scattering 
angles. 



C. The Case of an Axial Field B = B Zm 

For an axial magnetic field with B — Bz,(w) Z*, we have that ip s = and 9 S — i. The scattering phase 
functions are quite similar to the zero field case, except that now the Hanlc effect appears in factors in the 
functions C and D. Using Paper II, the phase functions are given by 

( hi = 1 + f-Bi \\ (1 -3cos 2 i) +sin 2 iC(ro,v?)] 

h=< h Q = |£a [sin*i- (l + cos 2 /) C{m,<p)] (CI) 

[ hu = jEi cos i D(m, ip). 

where 



where 



and 



C(n7, Lp) — cos 2 «2 cos 2<p — \ sin 2«2 sin 2<p (C2) 
D(zu,ip) = — cos 2 oli sin 2ip — i sin 2a^ cos 2ip, (C3) 



cos 2 a 2 = TT ^, (C4) 



cosa 2 smtt 2 = 757^ ■ (Co) 

1 "T Oq t 

Solutions for the vector Stokes flux is no longer analytic. With the Hanlc effect, T\j ^ except for B = 
or in the saturated limit. In the latter case of b 3> 1 everywhere, an analytic solution can be obtained, 
which is given by 



Jf c = (1 + cos 2i) C(m) (C6) 

3 
4 

Note that in this limit, the relative polarization becomes 



T S Q C = IE, sin 2 iJT(m). (C7) 



„s _ ™ / x-sc _ 3£i sin 2 i 

which is a constant across the profile and a function of viewing inclination only. This means that the polarized 
profile is flat-topped. 

Figure [10] displays profiles for the axial field case with E\ = 0.5 at a fixed value of sin 2 i = 0.4 but 
with different field values at the inner disk radius of w = 1, of logfoo = —2, —1, 0, +2, +4 to achieve a large 



1G 



dynamic range in Hanle ratios throughout the disk. As bo increases, the location where B = i?Han moves 
outward to w = bo. The upper left panel in Figure [8] shows normalized profiles of J-"f c ; lower panels show 
the fractional polarizations q s u and u s v ; and upper right shows the q s u -u s v shapes. The color sequencing is the 
same as in Figure [5] In large part the effect of an axial held is to rotate and foreshorten the Q-U segments 
relative to the zero field case. 



D. The Case of a Toroidal Field B = B v 

For a toroidal magnetic field with B = B v <p, the spherical geometry for the scattering problem is moder- 
ately complex. It is not possible to write down simple complete expressions for the phase scattering functions. 
But as demonstrated in Paper II, there are still some special cases that are analytic. For example in the 
saturated limit, Fff = 0, and the / and Q fluxes become 

(1 + \Bk) J*°(m) - §£i sin 2 i h FJ°(m - 1) (Dl) 
-§£i [j^°(m)-(l + cos 2 z)i J^°(m-l)]. (D2) 

For a disk viewed edge-on, solutions for the Stokes fluxes cannot be derived analytically; however, the 
scattering functions simplify to 

!hi = 1 + §-Ei (1-3 sin 2 <p) + |£i cos 2 ip cos 2 a 2 , 
h Q = -|£a cos 2 ip + §Ei (1 + sin 2 tp) cos 2 a 2 , (D3) 
hjj — sintp sin2a2- 

Examples of scattering and polarized profiles are shown in Figure [TT] as well as Q-U plots across the polarized 
profiles. In the previous section, sin 2 i = 0.4 was used for an axial field to give a significant signal in u%. For 
this toroidal field case, an edge-on disk with sin 2 i = 1.0 was used. The style of this figure is the same as 
Figure [TU] The a major distinctive in relation to a disk with an axial field is that a toroidal field leads to u s v 
profiles that are antisymmetric instead of symmetric. 
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Fig. 10. — Similar to Fig. [51 here for the case of an axial field in the point star approximation. 
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Fig. 11. — Similar to Fig. [TU1 here for the case of a toroidal field in the point star approximation. 
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